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r~' . Abstract 
t — . 

, The basic concepts and equations of Newtonian Cosmology are presented in the 

CN I form necessary for the derivation and analysis of the averaged Navier-Stokes-Poisson 

^ ■ equations. A particular attention is paid to the physical and cosmological hypothe- 

O . ses about the structure of Newtonian universes. The system of the Navier-Stokes- 

I Poisson equations governing the cosmological dynamics of Newtonian universes is 

O^' presented and discussed. A reformulation of the Navier-Stokes-Poisson equations in 

^ . terms of the fluid kinematic quantities is given and the structure of this system of 

equations is analyzed. 

> ■ 
^> . 

H ■ 1 Introduction 

In order to derive a system of the averaged Navier-Stokes-Poisson equations and to make 
use of the system for analysis of inhomogeneous Newtonian universes, it is significant to 
present the basic concepts and equations of Newtonian Cosmology in the form necessary 
to approach this problem. In this presentation a particular attention will be paid to the 
physical and cosmological hypotheses about the structure of Newtonian universes which 
are spatially infinite configurations of the self-gravitating compressible Newtonian cosmo- 
logical fluid without viscosity and heat transfer. The system of the Navier-Stokes-Poisson 
equations governing the cosmological dynamics of Newtonian universes is presented and 
discussed. A reformulation of the Navier-Stokes-Poisson equations in terms of the fluid 
kinematic quantities is given and the structure of this system of equations is analyzed. 

The structure of this paper is as follows. Chapter 2 describes the full set of the Navier- 
Stokes-Poisson equations for the self-gravitating compressible Newtonian fluid without vis- 
cosity and heat transfer under the gravitational field of its own Newtonian gravitational 
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potential. The definition of the Newtonian-space and analysis of the coordinate freedom 
allowed by the Kinematical group of the Navier-Stokes-Poisson equations is carried out 
in Chapter 3. It is pointed out, in particular, that the Newtonian gravitational potential 
is not uniquely defined for spatially unbounded fluid configurations. Chapter 4 is devoted 
to the discussion of the initial boundary value problem for the Navier-Stokes-Poisson 
equations which is not well-posed and the Heckmann-Schiicking boundary conditions for 
homogeneous and isotropic fluid configurations are formulated. The main hypotheses of 
Newtonian cosmology are formulated in Chapter 5. In Chapter 6 the Newtonian cosmo- 
logical principle is formulated and the class of the homogeneous and isotropic Newtonian 
universes is studied. A definition of inhomogeneous and/or anisotropic Newtonian uni- 
verses is given in Chapter 7. The nonexistence of the static homogeneous and isotropic 
Newtonian universe with the vanishing Newtonian cosmological constant is proved in 
Chapter 8. Chapter 9 makes a brief comparison of Newtonian and general relativistic 
cosmologies and gives evidence in favor of the remarkable physical similarity of the both 
cosmological pictures. The full system of the Navier-Stokes-Poisson equations in terms of 
kinematic quantities is derived in Chapter 10. The next Chapter gives a detailed discus- 
sion of this system of equations. It is pointed out, in particular, that in order to become 
closed it requires an additional evolution equation for the tidal force tensor which is con- 
structed from the gravitational potential and is an analogue of the space-time curvature 
tensor of general relativity. It means that the cosmological fluid configurations governed 
by the Navier-Stokes-Poisson equations in terms of kinematic quantities do evolve due to 
a generalized Newtonian gravitational force. In the last Chapter 12 the Raychaudhuri 
evolution equation of Newtonian cosmology is shown to be equivalent to the Friedmann 
equation for the cosmological scale factor for homogeneous and isotropic Newtonian uni- 
verses. 

The formulae and Sections from the paper I |T|| will be referred to as (I- XX) and I-X, 
correspondingly. 

Conventions and notations are as follows. All functions / = are defined on 3- 

dimensional Euclidean space in the Cartesian coordinates {a;*} with Latin space indices 

k, ... running from 1 to 3, and t is the time variable. The Levi-Civita symbol eijk is 
defined as £123 = +1 and e^^^ = +1 and 5* and 6ij are the Kronecker symbols. The 
symmetrization of indices of a tensor TJ^ = TjJ^x\ t) is denoted by round brackets, T^^^^ = 
I (Tj^ + Tfcj), and the antisymmetrization by square brackets, Tj*.^j = ^ (Tj^ — T^j). A 
partial derivative of Tj^^. with respect to a spatial coordinate or time t is denoted either 
by comma, or by the standard calculus notation, Tlj^ ^ = dT^j^/dt and Tj;. ^ = dT'!j^/dx^ , 
the fluid velocity is = u\x^,t), and the material (total) derivative is Tj^, = dT^i^/dt = 
dTji^/dt + u^dTji^/dx^ = T^j^ ^. + u^T'^^^ i. The Newton gravitational constant, the velocity 
of light and the Newtonian cosmological constant are G, c and A, respectively. 
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2 The Navier-Stokes-Poisson Equations 



Let us consider the self-gravitating compressible Newtonian fluid without viscosity and 
heat transfer under the gravitational fleld of its own Newtonian gravitational potential 
= 0(a;*, t). The motion of the fluid is determined by the fluid velocity field = u^{x\t) 
satisfying the compressibility condition (1-18), 

everywhere and by the fluid density fleld p — p{x\t) positive everywhere (1-29), 

p{x\t)>0. (2) 

The symmetric fluid stress tensor T*-' = T^^{x^,t) (1-40) satisfying the Boltzmann postu- 
late (1-75), T*-' = is assumed to be the fluid stress tensor of the perfect fluid, that 
is, the Newtonian fluid (1-83) with the vanishing viscosity coefficients, A(x\i) = and 
p{x\t) — 0, the vanishing heat fiux vector q^{x^,t) — and the constant temperature 
T = const, (1-79) such that 

rpij ^ _^ijp_ (3) 

An equation of state relating the fluid pressure p — p{x^,t) and the fluid density p — 
p{x\t) is assumed to be deflned (1-81) as 

The equation of state 

p = p{p) or p = p{p). (4) 



The external force = F^{x^,t) for the self-gravitating fluid is deflned to be due to 
the Newtonian gravitational potential 4> = 4>{x\t), 

F' = (5) 

where gi = gi{x^,t) is the Newtonian gravitational acceleration vector of a fluid particle 
moving with the velocity u^{x^ ,t). The acceleration vector has the following property. 

The identity for the Newtonian gravitational acceleration g^j, 

e'''g,,k = 0, (6) 



which means that the Newtonian gravitational field is always locally represented by its 
scalar potential (t>{x^,t) 

gi = (t>,i- (7) 

The Newtonian gravitational acceleration vector gi{x\t) and the Newtonian gravitational 
potential 4>{x\t) satisfy the Poisson equation with the Newtonian cosmological constant 
A 
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The Poisson equation for the Newtonian gravitational potential 

6''^gij = A-nGp - A, or 5*^0,^ = Ar^Gp - A. 



The conservation of fluid mass (1-30) brings the equation of continuity (1-31) for the fluid 
density p(x*, t) 

The equation of continuity for the fluid density p 



P,t+(pw*), = 0, (9) 



and the conservation of linear momentum (1-37) brings the Cauchy equation of motion 
(1-41) which for the fluid stress tensor (^) and the external force (^ becomes[| 

The Navier-Stokes equation of motion for the fluid velocity u' 

Ui^t + Uiju^ = -gj - ^Pj or + Uiju^ = -4>,i - -^P,j- (10) 

Another form of the Navier-Stokes equation (^) can be written by using the total fluid 
acceleration vector = A^{x\t) (1-45) for a fluid particle, 

dui 

Ai = — + gi = ai + gi. (11) 

The total fluid acceleration vector yl*(x-', t) can be defined (1-46) by the Navier-Stokes 
equation of motion (|T0|) as 

The total fluid acceleration Ai as determined by the Navier-Stokes equation 

Ai = --p,i. (12) 
P 



This form of the Navier-Stokes equation means from the physical point of view that 
the total acceleration Ai of a fluid particle due to the combined effects of the inertial and 
gravitational forces represented by the fluid particle's acceleration Oj and the Newtonian 
gravitational acceleration vector g^ is determined by the gradient of the fluid pressure. 



Due to the Navier-Stokes equation of motion (|T2D the total fluid particle acceleration { ^A\) 
is always directed away from a region with higher pressure towards a region with lower 
pressure. 



^ Though the equation of motion ( |lO[) is the Euler equation for the perfect fluid (1-80), it will be called 
here the Navier-Stokes equation of motion in accordance with Eq. (1-83) where the viscosity coefficients 
are assumed to vanish and the external force is due to the Newtonian gravitational potential (|^). 
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The system of seven first order partial differential equations (^), and ( ]TI1| ) has 

seven unknowns u\x^,t), p{x\t) and gi{x\t) with a given equation of state (H). It is 
equivalent to a system of five first and second order equations (||), and (|l3) for five 
unknowns u'^{x\t) , p{x^,t) and This system is of a mixed hyperbolic-elliptic 

type and in order to be solved in a physical setting of interest it should be supplemented 
by the corresponding initial and boundary conditions for the unknowns. As is has been 
discussed in Section 1-13 the problem of solving the equations of classical hydrodynamics 
is very difficult. The Poisson equation brings additional difficulties when one is look- 
ing for a solution for a self-gravitating fluid configuration satisfying the system of the 
Navier-Stokes-Poisson (||), (|) and ( 1To|) . A few results have been proved rigorously 
for compact Friedmann-like solutions and their finite perturbations for the polytropic 
equation of state, in particular, the linearization stability, a local-in-time existence of 
the solutions given a set of data not precisely equivalent to initial data, the occurrence of 
kinematical singularities (see and reference therein). Conditions for the equilibrium 

of the self-gravitating rotating systems for particular classes of the equations of state have 
been also studied (see p and references therein). 



3 The Newtonian Space-time and the Kinematical 
Group 

The Navier-Stokes-Poisson equations (^, (P) and (0) have been formulated in the frame- 
work of the field description of fluid motion when all objects characterizing the fluid and 
its properties are defined with respect to a position {x*} at a time t, see Section 1-5. Before 
analyzing and solving the equations one must analyze the mechanics of the fluid motion 
and reveal the freedom in the choice of possible Eulerian coordinate systems, which is 
available as far as the equations are formulated in this framework. 

As it has been pointed out in Section 1-9, the Cauchy equation of motion (1-41) is 
the Newton second law for a moving fluid particle. From the point of view of Newto- 
nian mechanics the dynamics of the whole fluid configuration under consideration can 
be represented in the field description by a sequence of its space configurations given as 
3-dimensional spaces represented by the full set of all fluid particle's positions {x*} which 
are occupied at each instant of time t. The time t at positions {x*} is assumed to be 
measured by a family of ideal clocks each of which indicates values of time regardless of 
the previous motion. This is the concept of absolute time. On the basis of this picture of 
the fluid motion one can naturally define the notion of a Newtonian space-time manifold 
as follows 0, ig. 



The Newtonian space-time A Newtonian space-time is defined as a 4- dimensional 
differentiable manifold of the following structure: 

(a) the real-valued absolute time t is determined up to a linear transformation t — * 

^ The concept of a Newtonian space-time is known to have led to a geometic formulation of Newtonian 
gravity in the framework of the pseudo-Riemanian geometry Q , , ||] . 
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at + b, t,a,b G R, t > 0, a > 0; 

(h ) the 3- dimensional spaces t = const are Euclidean spaces endowed with a coor- 
dinate system given by the Eulerian space coordinates {x*}; 

(c) for each value of time t there is only one 3-dimensional space E\ 

Let us consider now the Cauchy equation of motion (1-41) in the absence of all external 
and internal forces, 

pJl = as = 0, = 0. (13) 

which corresponds the state of free motion of each fluid particle. The following theorem 
takes place |0. 

Theorem 1 (The inertial coordinates and the Galilean group) If a fluid is in free 
motion (j7^ there always exists a global system of the Eulerian coordinates called 
the inertial coordinates which are defined up to the Galilean group of arbitrary linear 
transformation, 

x' Aix^ + BH + C\ t^at + b, (14) 
where Aj, B^ and are real and constant, and the matrix Aj is nonsingular. 



Proof. The group of Galilean transformations ([IJ) is found as the general solution of 
the Cauchy equation for free motion (0) for the Eulerian space coordinates {x*} along 
fluid particle paths (1-2), = x\t), 

IF = (1^' 

and the law of allowed transformations of the absolute time given above. QED 

Evidently the paths of fluid particles moving freely ( |13|) are straight lines in the New- 
tonian space-time of the fluid. 

In the presence of gravitation when a self-gravitating fluid moves in its own Newtonian 
potential = 0(x'^,t) in accordance with the Cauchy equation of motion (1-41) with the 
external force (^ and the fluid stress tensor T*-' = 0, 

^ = -S'j^ as F^ = -5%, T^ = 0. (16) 
dt dxJ ^' ^ ' 

one cannot determine a state of free motion in the sense ([T3|) and ([T5|) and deflne the class 



of inertial coordinates (|1J). It is natural, however, to deflne another class of privileged 
coordinates which corresponds to the state of free fall (|T6|) of fluid particles in the grav- 
itational fleld. It should be pointed out here that the deflnition of free fall assumes that 
the inertial and gravitational masses of any fluid particle are the same so that given the 
same gravitational fleld any particles move along the same paths if their initial positions 
and velocities were the same. The following theorem determines the corresponding class 
of coordinates and its transformation group 0, [jl0|-|12 . 
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Theorem 2 (The free fall coordinates and the Kinematical group) // a fluid is 
in free fall motion j^Tdj ) there always exists a global system of the Eulerian coordinates 
{x*,t} called the free fall coordinates which are defined up to the Kinematical group of 
arbitrary linear transformation, 

x' ^ Aix^ + D\t), t^at + b, (17) 

where A* is the real-valued nonsingular matrix Aj and D^{t) is a real-valued arbitrary 



function of timet. Under the change of coordinates ( |7^ j the gravitational potential (f){x'^ , t) 
undergoes the transformation 

(18) 

Proof. The group of kinematical transformations ( P!7D and transformations of the po- 
tential (|I8D are found as the general solution of the Cauchy equation for free fall motion 
([IGD for the Eulerian space coordinates {x*} of a fluid particle path (??), = a;*(t), 



-6^^^, (19) 



dt'^ dx^ ' 

and the law of allowed transformations of the absolute time given above. QED 



It should be pointed out here that the form of the Cauchy equation of motion (|T6|) 
remains the same in any free fall coordinates (0) and ([TsD . 

Thus, in the presence of gravitation the coordinates of a position {x*} of a free falling 
fluid particle are deflned up to the time-independent rotation x* — >■ A^^x^ at a given instant 
of time t, and up to the time-dependent translation x* — >■ x* + D^{t) possibly different for 
different moments of time t. The Newtonian potential is not a coordinate-free function, 
but rather deflned in dependence of a particular choice of the Eulerian coordinates {x\ t} 
and it changes due to (p!8[ ) upon their change. There is, however, a particular class of the 
self-gravitating fluid conflgurations called the isolated fluid conflgurations when a fluid 
occupies a bounded compact region with a unique common Newtonian potential ^ , ]11 
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Corollary 1 (The isolated fluid conflgurations) The Newtonian gravitational poten- 
tial (f) ^TBJ for the isolated fluid configurations satisfying the global boundary condition at 
spatial infinity 

0(x^^) ^ as ((5ijX^x^y^^ ^ oo (20) 

is uniquely defined, 

0^0- (21) 
Then the Kinematical group ( jl^ j reduces to the Galilean group (|7^j, 



D^{t) = BH + C\ (22) 
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For the case of isolated fluid configurations ( PPD one is able to distinguish in an invari- 
ant manner the inertial and gravitational parts of the total acceleration (pJ]). Indeed, 



the left-hand side of the Cauchy equation of motion for the free fall (|T^) which is an iner- 
tial acceleration a* and its right-hand side which is a Newtonian gravitational acceleration 
g'^ are kept invariant independently under the Kinematical group (O) with (22). 



4 The Initial Boundary Value Problem 

In general case of the noncompact fluid configurations when the self-gravitating fluid 
and therefore the gravitational field are extended all over a noncompact 3-dimensional 
Euclidean space such a boundary condition at spatial infinity (^Of ) cannot be imposed. 
Such cosmological fluid configurations are considered in modelling Newtonian universes 
represented by Newtonian space-times filled with the matter in the form of self-gravitating 
cosmological fluid distributions over in noncompact 3-dimensional Euclidean spaces 
evolving during an interval of time possibly infinite. In those cases one is forced to make 
a choice either to abandon the concept of the inertial Eulerian coordinates (|I^) or to 
announce the free fall coordinate systems ([T7|) inertial. The choice is usually made in 



favour of the latter [0], [|lTl-[|13[. From the physical point of view, such a choice means 
that all observers belonging the same 3-dimensional Euclidean space E^ a.t t = const are 
inertial though such inertial observers may accelerate relatively to each other when they 
move with the fluid due its gravitation. 

The initial boundary value problem for the system of the Navier-Stokes-Poisson equa- 
tions (H), (^), (|) and ( p!OD is known not to be well-posed. A 3-dimensional Euclidean 
space t = const is a characteristic surface of the system and it is not permitted to to 
set initial data on such a surface [|T^. As a result, the potential is known up to an arbi- 
trary harmonic function = ilj{x^,t) satisfying a boundary value problem for the Laplace 
equation = since it is not defined uniquely (|I^). Now with the partial time 

derivatives of the fluid density dp/dt and the fluid velocity du^/dt being determined from 
the equation of continuity @ and the Navier-Stokes equation of motion ([l0|) , the time 
derivative of the potential d(j)/dt must satisfy a boundary value problem for the Poisson 
equation (|[) without A term, 

S-^^ = ^^G^. (23) 
dx'dx^ dt dt ^ ' 

which gives rise again to an arbitrary harmonic function ipi = ipii^^, t) satisfying a bound- 
ary value problem for the Laplace equation S^^ {ipi,t) ,ij = due to (|T8|). Such a process 
eventually leads to an infinite number of boundary value problems. 

A physical manifestation of the problem with the bad-posedness of the initial value 
problem for the Navier-Stokes-Poisson equations can be understood if one considers the 
simplest possible noncompact fluid configuration, that is, a self-gravitating fluid with 
constant density, p = p{x\t) = const. On the basis of physical arguments one expects 
a homogeneous distribution of the self-gravitating fluid throughout the all 3-dimensional 
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spaces t = const, since a uniform, constant density does not show any particular point in 
the fluid distribution as preferable. As there is no preferable point in the fluid distribution, 
it is expected to move as a whole and there should be a particular inertial frame where 
the fluid is at the state of rest at a time t = t^. A natural initial condition for the fluid 
velocity therefore is u\x^ ,to) = . The Navier-Stokes-Poisson equations (H), (||) and ( |T0D 
in this case read 

= '"^^ - ^' ^ = -^' P«^) = const, nV,to) = 0. (24) 



The Poisson equation in the system (p^) can be solved in the spherical coordinates (r, 6, (p) 
and the resulting solution the Navier-Stokes-Poisson equations 

uUx^t) = --UnGp- A)xUt-to), 0(xO = - (47rGp - A) r^, (25) 
3 6 

in contradiction with the initial assumptions for the classical Poisson equation (H) when 
the Newtonian cosmological constant vanishes, A = 0. Indeed, first of all, the solution 
(pSD appears to distinguish the origin of an arbitrarily chosen coordinate system from 

■ 1/2 

other points since the gravitational potential vanishes, 0(a;*) = 0, at r = {5ijx'^x^) = 0. 
Secondly, the fluid cannot be at rest at a 3-dimensional space t = to = const , but rather 
moves with a constant acceleration a*(x^) = —AnGrx^/S away from the point with r = 
where the gravitational potential 0(a;*) = 0. 

If now one assumes the condition AnGp = A, the situation improves, since in this case 
the solution (pSj) shows that u'^{x^,t) = and 0(a;*) = everywhere in agreement with 
the physical analysis made before approaching the equations. It should be emphasized 
here that the presence of the Newtonian cosmological constant A in the Poisson equation 
means a modification of the classical Poisson equation of Newtonian gravity where A = 0. 
This issue is discussed later in this Section. 

Therefore, any way to improve the situation with the initial boundary value problem 
for the Poisson equation and, as a result, with the Navier-Stokes-Poisson equations, to 
allow a proper analysis of self-gravitating fluid configurations defined on noncompact 3- 
dimensional spaces t = const, would assume a modification of the Poisson equation either 
in its structure, or in the structure of its boundary conditions. The latter also would result 
in an effective change in structure of the Poisson equation, as any boundary condition 
on the gravitational potential itself 0(x^,t) could not improve the situation as it is clear 
from the above analysis. 

The attempts of the first kind are being made since a long ago, for the unsatisfactory 
situation with the application of the classical Poisson equation to unbounded systems 
has been realized very soon after the discovery of this equation. In 1895 Seeliger and 
Neumann [jl5[, proposed the modified Poisson equation 

6'^^^-X<P = ^rrGp, x = const, (26) 

which assumes a finite range x~^^'^ of this gravitational force as compared with an infinite 
range of the Newtonian gravitational potential. This equation can be shown to have the 
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constant solution for the constant fluid density, p = p{x^,t) = const, 



uHx\ t) = 0, = —Gp = const, (27) 

X 

which in agreement with the expected homogeneous fluid distribution in the state of rest 
once it was at rest initially. However, one can show that this equation can not be obtained 
as a weak field limit of the Einstein equations of general relativity even in the presence 
of the cosmological constant A = Ac^ where A is the cosmological constant entering the 
Einstein equations. It is the the Poisson equation (|^) that follows as the weak field limit of 
the Einstein equations. This limiting procedure is singular in the sense that the Einstein 
equations for perfect fluids with an equation of state (^) have a well-posed initial boundary 
problem [|l7l, [|l^, but the Poisson equation has not. It has been assumed therefore that a 
modification of the Poisson equation must be dictated by a reconsideration of its derivation 
from the Einstein equations by taking into account the post-Newtonian approximation 
terms which provide the weak-field correction terms next to the Newtonian approximation. 
It has been shown []T9| that one can construct a generalization of the Navier-Stokes- Poisson 



equations which does possess a well-defined initial boundary problem, but those equations 
have very different structure and nonlinear in the gravitational potential itself and in 
its products with the fluid velocity. Another generalization of the Navier-Stokes- 



Poisson equations gives a system of equations very similar to the linearized Einstein 
equations which provide an equation for propagation of the gravitational potential, to 
permit gravitational waves and a finite speed of gravitational interaction as compared with 
the infinite speed propagating Newtonian gravitation with no gravitational radiation. 

Attempts to modify the boundary conditions for the gravitational potential have 
proved to be more fruitful because in this approach one can usually assume a condition 



on the basis of physical arguments |2T[ . First of all, it should be noted that using a simple 

'1/2 

generalization of the boundary condition (^) such as 0(^)|oo as {5ijX^x^) oo 
or a prescription for the value of lim^'^'^^fci as {6ijX^x^)^^'^ oo, does not work because 
the potential (p diverges at spatial infinity and ^^^(p^ki is determined by the limiting fluid 
density at infinity though the Poisson equation. Indeed, the structure of the Poisson 
equation (H) is known to determine the divergence of the acceleration, 5^^gij = S^^cp^ij, 



while the so-called tidal force tensor Ey = Eij(x ,t) defined as 

The tidal force tensor Eij for the Newtonian gravitational potential (p 

Eij = 4>^ij — -5ij5^''(j)^ku Eij = Eji, 5^^Eki = 0, (28) 



is left undetermined, the tensor being an analogue of the space-time curvature of general 
relativity. Therefore eight component of the trace-free symmetric tensor Eij are not de- 
termined by the Poisson equation and they are therefore pure kinematic quantities in the 
Newtonian gravity. The absence of the field equation to govern the tidal force tensor Eij 
(pSD may be considered as a reflection of the above problems. On the other hand, one can 
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use the freedom in Eij to get a suitable boundary conditions for particular classes of spa- 
tially noncompact self-gravitating fluid configurations. The suitable boundary conditions 
that allow one to handle spatially homogeneous fluid configurations which necessarily 
have the uniform fluid density have been formulated by Heckmann and Schiicking |Tl| , 



12 , 21 



The Heckmann-Schiicking boundary condition for the tidal force tensor E^j 

The following boundary condition is assumed to hold for spatially homogeneous and isotro- 
pic fluid configurations 

Eij{x^,t) ^ as [Sijx'x^Y''^ oo. (29) 



This boundary condition ( pof ) enables one to describe consistently the unbounded homoge- 
neous and isotropic distributions of the self-gravitating fluid. When taken as cosmological 
models, such cosmological fluid configurations are Newtonian analogues 0, [|lll]-[|13] of the 
Friedmann-Lemaitre- Robertson- Walker (FLRW) cosmological models of general relativity 
p2[-p6[. The conditions, however, can be shown to exclude the spatially homogeneous 



but anisotropic self-gravitating cosmological fluid configurations. Such Newtonian con- 
figurations corresponds to the Kantowski-Sachs and Bianchi spatially homogeneous but 
anisotropic cosmological models of general relativity |26|-|2^. A generalization of the 
boundary condition (^) to cover this broader class of fluid conflgurations is also known 

a Hi. 0- 



The generalized Heckmann-Schiicking boundary condition for the tidal force 
tensor Eij The following boundary condition is assumed to hold for spatially homogeneous 
anisotropic fluid configurations 



,t)^E{t) 



ij I oo 



as 



1/2 



OO. 



(30) 



Here E{t)ij\oo are arbitrary functions of time. The boundary conditions (p9|) and ( pO]) 
are considered as physically adequate for this class of cosmological fluid conflgurations 
(see for a discussion and further references). In both cases information immediately 



propagates in from inflnity to determine the local physical evolution of a moving fluid. 
In this approach no new equation for the propagation of the tidal force tensor Eij is 
explicitly necessary if one remains in the framework of Newtonian gravity, though it leads 
to impossibility to decide a priori which boundary condition flts a physical setting under 
consideration. If such a propagation equation for dEij/dt is deflned in a framework of 
the generalized Newtonian gravity, then a physical setting would dictate a particular 
boundary condition at inflnity, to determine nonlocally through this condition the local 
dynamics of a self-gravitating cosmological fluid. 
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5 The Newtonian Cosmology 



Cosmology is a theory of the evolution and structure of our Universe. In dependence on 
physical and mathematical hypotheses about the structure of the Universe space-time, 
the content and distribution of cosmological matter and the set of dynamical equations 
governing the Universe evolution one can formulate different theoretical frameworks for 
derivation and analysis of cosmological models to compare their results with observations 
and experimental data (see 0, |13|, p9|- for a discussion and references). The 



following hypotheses are assumed for Newtonian universes, the cosmological models of 
Newtonian Cosmology. 

The Newtonian cosmological space-time hypothesis Newtonian universes are as- 
sumed to have the geometry of Newtonian space-time. 

A Newtonian universe therefore always possesses an absolute time t. A state of a 
Newtonian universe at every instant of the absolute time t is given by a 3-dimensional 
Euclidean space and can be therefore described by the Eulerian coordinates {x^,t) 
(see Section 

The Newtonian cosmological matter hypothesis A Newtonian universe is assumed 
to be filled with a self-gravitating compressible Newtonian cosmological fluid without vis- 
cosity and heat transfer under the gravitational field of its own Newtonian gravitational 
potential. 

Though generally a cosmological fluid can be assumed Stokean, see Section 1-13, the 
Newtonian cosmological fluid is taken as a physically reasonable model for the matter in 
our Universe and due to its own gravitation as the only reason of the universe evolution. 
The next hypothesis determines the physical laws of fluid dynamics and its gravitation. 

The Newtonian cosmological dynamics hypothesis A Newtonian universe is as- 
sumed to be governed by the Navier-Stokes-Poisson equations ij^) and ^Tdj) sup- 
plemented by a set of initial and boundary conditions for a cosmological fluid configuration 
under study. 

Thus, a particular cosmological fluid conflguration determines a particular Newtonian 
universe. A particular cosmological fluid velocity fleld u^{x\t), the cosmological fluid 
density distribution p{x^,t) and the Newtonian gravitational potential fleld (^{x^ ,t) given 
by a solution to the system of the Navier-Stokes-Poisson equations (^, (H), (||) and ([1^ ) 
for a Newtonian cosmological model determine the evolution of the Newtonian universe 
as that of the cosmological fluid conflguration. The fluid velocity fleld u^{x\t) determines 
the equation of motion of cosmological fluid particles = x*(^-', t) and their trajectories 
(I-l)-(I-6) through the differential equation (I-IO), 

— = u\x^,t), x\0)=e, (31) 

to predict the state of the Newtonian universe for moments of the time interval of the 
cosmological evolution allowed by the initial and boundary conditions in dependence on 
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the initial positions x*(0) = ^' of the cosmological fluid particles. With a particular 
distribution of the cosmological fluid density p{x^,t), the motion of the fluid particles 
taken as matter constituents should be compared with the observed motion of the real 
cosmological matter to conclude about the physical adequacy of the cosmological model 
and interpret and/or predict the data. 



6 The Newtonian Cosmological Principle 

The cosmological matter in our Universe is distributed in a highly inhomogeneous manner 
on the scales characteristic of typical matter condensations such as stars, galaxies, cluster 
of galaxies, etc. However, for the largest scales the overall distributions of the matter 



structures is assumed to be homogeneous and isotropic |jT3[, ||2^- ||3T|]. This fundamental 
assumption about the large-scale structure of the real Universe is called the Cosmological 
Principle. Let us formulate this Principle in the framework of Newtonian Cosmology , 

0, m, 0. 

The Newtonian cosmological principle For any observer moving with a cosmological 
fluid particle the Newtonian universe appears to be the same at different points of space, 
at different instances and at different directions. Any other observer moving with another 
cosmological fluid particle at the same time observes the same appearance of the Newtonian 
universe. 

The Newtonian cosmological principle assumes the existence of a privileged family of 
observers who, at each instant of time, observe the Newtonian universe as homogeneous 
and isotropic. The distributions of the cosmological fluid over 3-dimensional Euclidean 
spaces at each instant of time are therefore homogeneous and isotropic. However, that 
does not mean that the cosmological fluid is static. In general, it does evolve in time since 
the structure of the Newtonian space-time permits nontrivial time shifts (|1^) between the 
3-dimensinal Euclidean spaces t = const. A fundamental theorem based on the result of 



McCrea and Milne ||35|, establishes the law of the evolution of a Newtonian universe 
satisfying the Newtonian cosmological principle 0, [|l^. Up to a different interpretation 
of one term, this equation has the same form as the Friedmann equation [^] governing 
the evolution of the FLRW universes which are homogeneous and isotropic cosmological 
models in the framework of general relativity [|r^, ||S^, P7|. 



Theorem 3 (The homogeneous and isotropic Newtonian universes) The evolu- 
tion of a Newtonian universe satisfying the Newtonian cosmological principle is governed 
by the Friedmann equation 

UJ =^^°-^+3^^ ^''^ 

where R = R{t) is the so-called scale factor of the cosmological fluid, k is the constant 
of integration which represents the total energy of fluid particles, po = p{to) is the initial 
value of the cosmological fluid density. 
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Proof. Let us consider a Newtonian universe filled with a cosmological fluid satisfying 
the Newtonian cosmological principle an governed by the Navier-Stokes-Poisson equations 
), (P) and (p!oD. An observer moving with a fluid particle and having set up a 



system of the Eulerian coordinates {x^,t) with the origin at the fluid particle measures 
the physical properties of a fluid particle at a position {x*^} as its the velocity 
the density p{x^,t) and the pressure p{x^,t). Another observer residing on another fluid 
particle with another system of the Eulerian coordinates {x'^,t) with the origin at that 
fluid particle measures the physical properties of the same fluid particle at the same 
position as the first observer, to obtain the velocity t), the density p(x'^, t) and the 

pressure p{x'^, t). Dashes over the density and the pressure have been omitted since these 
fluid properties are defined independently of any observer. The Newtonian cosmological 
principle demands now that u"'{x'^,t), p{x'^,t) and p{x"',t) should be the same functions 
of x"' and t as u^{x'',t), p{x^ ^t) and p{x'^,t) are functions of x^ and t, 

t) = t), p(x'^ t) = p(x^ t), t) = p{x\ t). (33) 



The second and third conditions ( P^ ) mean that the density and pressure should be 
independent of a position of a fluid particle, 

p = p{t), p = p{t), (34) 

and they are functions of the absolute time t only. If the relative velocity of the fluid 
particle associated with the first observer with respect to the fluid particle associated 
with the second observer is u'^{x'^ — x^,t), then, by taking into account the first condition 
(|33D , the velocity u'\x'^,t) of the fluid particle at the position {x'^} can be represented 
as 

u\x'\t) = u\x\t) + u\x'^ - x\t) or u\x'^ -x^,t) = u\x'\t)-u\x^,t). (35) 

The system of functional equations ( ^5|) has a solution 

u\x\t) = Hi{t)xK (36) 

Now by the Newtonian cosmological principle the motion of the cosmological fluid is 
isotropic, that is, it must not depend on any direction, -f^j(t) = H{t)6j, and the cosmo- 
logical fluid velocity (^) takes the form 

u'{t) = H(t)x\ (37) 

If {.^*} is an initial position of a cosmological particle at t = to? 2;* (to) = ■CS the solution 
to this initial value problem gives the equations of motion of a cosmological fluid particle, 
see Eq. (0), in the following form 

1 flT-^ 

x\t) = R{t)C where -— = H{t), R{to) = 1. (38) 
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This shows that the only motions of the cosmological fluid compatible with the Newtonian 
cosmological principle, that is, with a homogeneous and isotropic distribution of the fluid, 
are those of the uniform expansion or contraction (|38D determined by one time-dependent 
scalar scale factor R(t) which is always positive, R(t) ^ 0. This competes the analysis of 
the cosmological fluid kinematics due to Newtonian cosmological principle. 

To study the dynamics of the cosmological fluid satisfying the Newtonian cosmological 
principle, one should analyze the system of the Navier-Stokes-Poisson equations (^), 
(0) and dTDI ) for the fluid moving in accordance with Eqs. (^Tj) and (pS]). The equation of 
continuity (|) takes the form 

1^ + 3H{t)p{t) = ^ + 3H{t)p{t) = 0, (39) 
which has a general solution 

m - (40) 

The solution (|40|) admits an obvious interpretation, since if all linear dimensions in an 
evolving homogeneous and isotropic cosmological fluid are scaled up by the factor R(t), 
then all volumes are scaled up by the factor R^{t) with the correspondent increase or 
decrease in the cosmological fluid density. 

The Poisson equation for the gravitational potential = (j){x^ ,t) of a homogeneous 
and isotropic distribution of the cosmological fluid with the density p = p{t) and with the 
Heckmann-Schiicking boundary condition (pQ]) reads 



a2 

S' \ . 0(x^t) = 47cGp(t) - A, E,Ax\t) ^ as idi^x'x^Y'^ oo. (41) 
ox'^ox^ 

This equation has a solution 

4>{x'',t) = '^'nGp{t)6ijx'-x^ - l-ASijX^x\ and (p^ij = ^TTGp(t)Sij - ^ASij. (42) 



The second relation in (|42|) shows that the Heckmann-Schii eking boundary condition ( pOj) 
is satisfied by the solution and, in fact, the tidal force tensor Eij (|28|) vanishes everywhere 

Eij = 0. (43) 
The physical meaning of the solution ( ^2]) is that the gravitational potential (j){x^,t) 



is spherically symmetric and constant over surfaces SijX^x^ = const for each moment 
t = const. Now the Navier-Stokes equation (|T0[) takes the form 

du' ,du^ 1 ^ii dp ^ii d(j) f dH 4 • 1. 

+ + = {-Jr + H']x' = -^^Gp{t)x^ + -Ax\ (44) 



dt dx^ p dx^ dx^ \ dt J 3 3 

By expressing the function H{t) and the fluid density p{t) in (|4^) through the scale factor 
R(t) by (^) and (^0]) one gets the equation for R(t), 

+ l^Gpito) - ^AR' = 0. (45) 
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The first integral of this equation (^) can be easily found by integration, to bring the 
Friedmann equation (|32D for the scale factor R{t). Here —k is the integration constant 
which has meaning of the total energy of a cosmological fluid particle with a coordinate 
R{t) moving along its path (|38|). Indeed the integration of (^) gives 

where K{t) and V{t) are the kinetic and potential energies of the fluid particle. QED 



7 Inhomogeneous and Anisotropic Newtonian Uni- 
verses 

On the basis of the notion of homogeneous and isotropic Newtonian universe one can 
define a notion of an inhomogeneous and/or anisotropic Newtonian universe. 

Inhomogeneous and/or anisotropic Newtonian universes A Newtonian universe 
is inhomogeneous and/or anisotropic if it does not satisfy the Newtonian cosmological 
principle, for, at least, an interval of the cosmological evolution time. 

As it has been pointed out the Newtonian cosmological principle assumes the existence 
of a privileged family of observers who, at each instant of time, observe the Newtonian 
universe as homogeneous and isotropic. If a Newtonian universe is inhomogeneous and/or 
anisotropic, a free falling observer comoving with a fluid particle will be measuring de- 
viations in the velocities, density an pressure for a fluid particle under observation as it 
moves, and different values of these quantities for different fluid particles moving around 
the observer. From the geometrical point of view that means that distributions of the 
evolving cosmological fluid over 3-dimensional Euclidean spaces at each instant of time are 
different and therefore the fluid is inhomogeneous and/or anisotropic during its evolution. 



8 The Static Newtonian Universe 



The equation ( P2| ) may be integrated for different values of A and k to bring a number 
of Newtonian cosmological models similar to the homogeneous and isotropic cosmological 
models of general relativity 0, |]l3l, |3^. One of the most important results of the analysis 



is that if the Newtonian cosmological constant vanishes, A = 0, then no static Newtonian 
universes, that is, either spatially compact or noncompact self-gravitating distributions 
of cosmological fluid for which scale factor R(t) ( P?]) and (^) is constant, 

R{t) = const x'{t) = constf , u\t) = 0, (47) 



can exist [0, [|T3], ^^^^ Section ^ Therefore, Newtonian cosmology, as well as 

general relativistic cosmology, also predicts that our Universe is evolving in time and there 
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was a time when the evolution began. If the Newtonian cosmological constant does not 
vanish, A 7^ 0, then there exists a critical value of the Newtonian cosmological constant 
A = Ac for which the Newtonian universe is static and the cosmological fluid configuration 
in this case neither expands, nor contracts. This Newtonian universe is an analogue of 
Einstein's static cosmological model of general relativity. 

Corollary 2 (No static Newtonian universe with A = 0) No static Newtonian uni- 
verse exists if the cosmological constant vanishes, A = 0. A static Newtonian universe 
exists for k > if the Newtonian cosmological constant A and the scale factor R{t) have 
the critical values, Ac = k^/lGn'^G'^pl and Rc = AnGp^/k. 

Proof. Let us assume that for a homogeneous and isotropic Newtonian universe the 



Newtonian cosmological constant vanishes, A = 0. When the equation of motion ( ^Sj) for 
the scale factor R{t) of a fluid particle reads 

R''^ + ^^G'p(to) = 0. (48) 
Under the condition of staticity p7|) Eq. (|48| ) gives 



p(to) = (49) 

that in contradiction with the definition of a Newtonian cosmology as a configuration of 
a self-gravitating cosmological fluid with everywhere nonvanishing positive fluid density 
(^) and the positivity of the mass (1-30) of any fluid region S(t). Therefore, such a static 
fluid conflguration is impossible. 

Let us now consider the Friedmann equation ( p2D for the scale factor R{t). For a 
static Newtonian universe (^) the equation takes the form 

SttG ,1.9 

-k + -AR' = 0. (50) 

For A; > when the potential energy of a fluid particle is larger than its kinetic energy 
(|46D , the critical value of the cosmological constant, A^ = fc^/lGvr^G^Po' ensures that the 
function in the right-hand side of Eq. ( |5(]| ) always positive for all R{t) ^ except it 
vanishes at a double root at R{t) = const = 4:7rGpo/k = Rc- One can show that for the 
cases k = and k < this function never vanishes for R{t) ^ 0. QED 



9 The Newtonian Cosmology versus the General Rel- 
ativistic Cosmology 

The general relativistic Friedmann equation for the scale factor RQn{t) has the same form 
( P^ with three differences: 

[a] the pressure is absent in the case of general relativity; 
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[b] the Newtonian scale factor R (t) has no physical dimension while general relativistic 
scale factor -Rgr(^) has that of time; 

[c] the constant k which is an integration constant representing the total energy 
—k (^6|) in Newtonian cosmology, in general relativity represents the curvature of 3- 



dimensional space orthogonal to the world lines of cosmological fluid particles, k having 
one of the value +1, or —1 for three different types of 3-dimensional space geometry. 

In all other respects the identity of the Newtonian and general-relativistic Friedmann 
equation ( P^ ) is complete |]TB[, [^, [^. There are three physical consequences of this 



remarkable similarity of fundamental significance for the construction of a realistic homo- 
geneous and isotropic cosmological model of our Universe in the frameworks of Newtonian 
and general relativistic cosmologies. 

{a} The history of an homogeneous and isotropic Universe as described in the frame- 
works of Newtonian and general relativistic cosmologies is identical for both theories since 
the Friedmann equation (|32D for the variation of the scale factor dR/dt has the same struc- 
ture and there is the same set of cosmological models in both cosmological settings. 

{6} The difference between the Newtonian and general relativistic cosmologies is gov- 
erned by the [pressure] / [density] and [gravitational potential] /[rest mass] ratios. Therefore 
for the current state of our Universe where both ratios are always small, general relativity 
cannot offer anything radically new. When the propagation of light over the distances 
comparable with the radius of the curvature of Universe is to be taken into considera- 
tion, general relativistic cosmology will be more physically adequate compared with a 
Newtonian cosmology consideration. However, one cannot expect any radical differences 
in both approaches for the scales characteristic, for example, for galaxies in analysis of 
galaxy formation, galaxy structure, etc. 

{c} The ability of general relativity to deal with the cosmological settings when the 
[pressure] / [density] and [gravitational potential] /[rest mass] ratios are not appreciable is 
not much greater as compared with the Newtonian cosmology, though the results more eas- 
ily and directly interpreted in the framework of general relativity. On the other hand, the 
mathematical difficulties one encounters while solving cosmological problems in general- 
relativistic setting frequently obscure the physical interpretation of obtained results. 

10 The Navier-Stokes-Poisson Equations in Kinema- 
tic Quantities 

For analysis of the dynamics of inhomogeneous Newtonian universes, that is, the evolution 



of self-gravitating cosmological fluids, it is very useful to reformulate the system of the 
Navier-Stokes-Poisson equations (||), (H), © and (|l^) in terms of the kinematic quantities, 
see Section I-IO. 
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10.1 The Evolution Equations 



The first issue in the reformulation of the system is to derive the evolution equations 
for the fluid expansion scalar = 9{x\t) (1-59), the fluid shear tensor cxjj = aij{x'',t) 
(1-58), and the fluid vorticity tensor Uij = ujij{x'',t) (1-49) or the fluid vorticity vector 
tu* = C(j*(x-', t) (1-50) which must replace the Navier-Stokes equation of motion (|10| ) 
describing the evolution of the fluid velocity vector = u'^{x\t) ^3 . 



Theorem 4 (The evolution equations for the kinematic quantities) The Navier- 
Stokes equation ( ff^ j for the velocity vector and the Poisson equation ^ lead to the 
following system of evolution equations for the fluid expansion scalar 6{x'^,t), the fluid 
shear tensor aij{x^, t) and fluid vorticity vector uj^{x^ , t): 

The Raychaudhuri evolution equation for the expansion scalar 9 

df) 1 

^ + -e^ + 2(a^ - uj^) + ATiGp - A - 6'' A, = 0, (51) 
dt 3 

The propagation equation for the shear tensor aij 



The propagation equation for the vorticity vector cj* 

duj'^ 2 1 

— + -Ou' - S'^a^kOJ^ - -e'^^A, k = 0. (53) 
dt 3 ^"2 ^' ^ ^ 



Proof. To derive the evolution equations for the kinematic quantities, one needs to 
consider the tensor Uij which describes the spatial change of the fluid particle's velocity 
by taking a spatial derivative of the velocity vector Ui = SijuK Due to the Cauchy 
decomposition theorem, see Section I-IO, the tensor Ui^j can be always represented (1-46) 
in terms of the fluid expansion scalar 6 (1-59), the fluid shear tensor aij (1-58) and fluid 
vorticity tensor ojij (1-49) as 

The kinematic decomposition of the tensor Ujj 

Mij = cFij + -5ijd + Uij. (54) 



By applying two identities for the commutation of partial derivatives of Uij with respect 
to t and 
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The first identity for the tensor Uij 

'^i,jt = Ui^tj, (55) 

The second identity for the tensor Uij 

Ui,jk = Ui^kj, (56) 



to a spatial derivative of the Navier-Stokes equation (0) written with using the total 
acceleration vector (ITT]) for a fluid particle 



Ui^tj + ^ Ui,kjUi + S Ui^kUij + gij - Aij = 0, (57) 

one gets the equation 

Uijt + S'^^UijkUi + +9i,j - Aij = ^''^^i^kUij + - = 0. (58) 

Making use of the definition of the tidal force tensor Eij (|28|) , the Poisson equation (H) and 
substituting the decomposition rule (|^) into the spatial derivative of the Navier-Stokes 
equation (^) with taking consequently its trace, the symmetric trace-free part and the 
antisymmetric part using the definition for the vorticity vector (1-50) brings the the 
evolution equations for the fiuid expansion scalar 9 (0), the fiuid shear tensor aij (^), 
and fiuid vorticity vector cu* (^ ). QED 

It should be pointed out here that Eq. (|58| ) is the propagation equation for the tensor 
Uij along the fiuid particle trajectories. In absence of the fiuid pressure when the total 
fiuid acceleration vector = (|T2|), if one considers two neighboring free falling fiuid 



particles moving along their trajectories the equation (|58|) describes a relative change in 
the positions of both particles due to action of gravitation through the second spatial 
derivative of potential (p^ij. This is the Newtonian analogue of the geodesic deviation 
equation in general relativity [38| and the quantity (p^ij or the tidal force tensor Eij (^), 



is the Newtonian analogue of the space-time curvature tensor of General Relativity ||^, 



10.2 The Constraint Equations 



The identities (pq ) and (^) for the tensor Newtonian analogues of the so-called 

Ricci identities in general relativity, which relate the second derivatives of an arbitrary 
vector field with the space-time curvature tensor [Q. 

Due to the second identity (^) for the tensor Uij there is also another set of equations 
which put constraints on the kinematic quantities [3^. 
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Theorem 5 (The constraints on the kinematic quantities) The second identity 
for the tensor Ui^j ( [5^ j puts the following constraints on the fluid expansion scalar 6 {x'' ,t) , 
the fluid shear tensor aij{x'', t) and fluid vorticity vector u!''{x\ t): 



The first constraint equation 



2 

5^^{aij^k - uJij,k) - -Oj = 0, (59) 



The second constraint equation 

= 0, (60) 



The third constraint equation 

6^^'ia,k,i + u;,k,i)s"'^'' = 0. (61) 



Proof. The first constraint (|59|) follows from taking a trace of the second identity ( ^q) 
for the tensor Uij with respect to indices i and j. The second constraint ( |60D follows from 
the fact that the total antisymmetrization of (|56D in its indices gives 



U[i,jk] = (62) 
which provides (pDj) after total contraction with the Levi-Civita symbol e*-'*^. The third 



constraint (|6TD follows after the contraction of the second identity for the tensor Uij ( ^6]) 
with e"^^^ and symmetrization of the two remaining free indices 

^^(n^m)Jk ^^^^^^ _ ^ q_ ^gg) 

QED 

It should be noted here that the trace of the constraint (^ ) vanishes identically due 
to the constraint (j^. 



10.3 The Integrability Conditions 

There are three more constraints on the kinematic quantities following from the first and 
the second identities (^) and (|56D for the tensor Uij as their integrability conditions. 

Theorem 6 (The intergrabihty conditions for the kinematic quantities) There 
is the following set of the integrability conditions of the first identity ([5^ and the sec- 
ond identity ( |5^ j for the tensor Uij . No more integrability conditions exist. 
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The first integrabifity condition 



„i SirG 

Kj = ^P.^ (64) 



The second integrability condition 

eI:/^'^ = 0, (65) 

The third integrabihty condition 

4a + li^fl = 0- (66) 

Proof. The derivation of the integrability conditions is straightforward. The first and 
the second intergrability conditions follow from taking a spatial derivative of the first 
identity (^) and the requirement Ui,t[jk] = which brings an integrability condition 

Wjjtfc — Ui,ktj ~ "^i.tb'fc] = ~^ Ui^jtk — Ui^ktj = 0. (67) 

Taking the trace of (^) with respect to indices i and j with using the evolution equations 



and the Poisson equation (§) gives the first integrability condition (|6^), while the antisym- 
metrization of ( pTf ) with respect to indices i and j gives the second integrability condition 
(|65D. No more constraints on the kinematic quantities can be found from the integrability 
condition (|67D . The third intergrability conditions follows from taking a spatial derivative 
of the second identity (|56|) and the requirement Uij[ki] = which brings an integrability 
condition 

Ui^kji — Uijjk — Uij[ki] =0 Ui^kji — Ui^ijk = 0. (68) 

Taking the antisymmetrization of ( |68D with respect to indices i and j gives the third 
integrability condition (|66D. No more constraints on the kinematic quantities can be 
found from the integrability condition (|68D. 

It is easy to show that calculations of the integrability conditions for ( |67D and ( |68D 
does not bring any new constraints since they are satisfied identically. QED 



It should be noted here that the trace of the integrability condition ( p3| ) vanishes 
identically because the tensor Eij is symmetric. 

The integrability conditions (0), (^ ) and ( pB| ) are analogues of the Bianchi identities 
for the space-time curvature tensor of general relativity |Q that can be seen from the 



structure of the integrability conditions ( pTl) and (^) which are just a cyclic combination 
for the third derivatives of the fluid velocity m*. Therefore taking further integrability 
conditions results in expressions identically vanishing in a 3-dimensional space for they 
involve the antisymmetrization of more than three indices. 
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11 The Structure of the System of Equations 



While the Poisson equation (P) keeps the same form in this reformulation of the Navier- 
Stokes-Poisson equations in terms of the kinematic quantities, the equation of continuity 
(^) takes the following form. 

The equation of continuity as the evolution equation for the fluid density p 

^ + P0 = 0. (69) 



Thus, the system of the Navier-Stokes-Poisson equations for the fluid expansion scalar 
6{x'^,t), the fluid shear tensor (Tij{x^^t) and fluid vorticity vector uj'^{x\t) includes the 
evolution equations (^T])-(^3p, the constraints (pUD-(|BTD, the integrability conditions 



(|66D , the equation of continuity (^), the Poisson equation (H) and (H), to replace the 
Navier-Stokes-Poisson equations (^, (^, (|^), (|^) and (PUf). Upon having determined the 
kinematic quantities by solving the above system, one must solve the nonhomogeneous 
equations for the fluid velocity u^{x^,t) together with the identities (0) and (^). 
Finally, to determine the equation of motion of fluid particles = x^{^\t), or C,-' = 
$,^{x\t), (I-l)-(I-4) from the fluid velocity u\x\t) the initial value problem ( Pl) has to 
be solved. 

This system of equations for the kinematic quantities and the Navier-Stokes-Poisson 
equations are usually considered to be equivalent [Q. However, this has not been shown 
as yet, strictly speaking. As it has been noted above, the Navier-Stokes-Poisson equations 
is a system of a mixed hyperbolic-elliptic type for seven first order nonlinear partial 
differential equations for seven unknowns p, m* and Qi, or it can be taken as a system of 
five first and second order equations (|p, and (|TU]) for five unknowns p, and cf) to 
be supplemented by the corresponding initial and boundary conditions for the unknowns 
and an equation of state It allows application of some techniques for analyzing the 
questions of existence and uniqueness, at least, for a class of boundary and/or initial 
conditions, see Section |[ Once reformulated in terms of kinematic quantities, the system 
of nonlinear equations for a self-gravitating fiuid changes its structure and it is does not 
fit into the standard classification of the partial differential equations. Four main issues 
appear to be important in consideration of the Navier-Stokes-Poisson equations written 
in terms of kinematic quantities: 

(A) The system is not closed since the evolution equation for the shear tensor ( |5^ ) con- 
tains the tidal force tensor Eij (^8]) which does not have any evolution equation for itself. 
This is a refiection of the noncloseness of the Navier-Stokes-Poisson equations (^), (^), 
(H) and ( |I0| ) with respect to boundary conditions and its initial boundary value problem, 
see Section ^. Upon reformulation in terms of kinematic quantities this problem becomes, 
however, critical for the Navier-Stokes-Poisson equations. Indeed, the system (||), (y) and 
(|10D can be supplemented by some boundary conditions, for example, the Heckmann- 
Schiicking boundary conditions ( PPj) or (POD, to allow its solution without adding any 
other equations and an explicit change in the structure of the system of equations (H), (0 
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), (g> and (|TU|). On the contrary, the system (g), (g), (|g), (g)-(g), (0)-® 
and (|69D must be exphcitly supplemented by an evolution equation for the tidal force 
tensor Eij (|^) in order to allow its solution. Thus, the Navier-Stokes-Poisson equations 
in terms of kinematic quantities must admit a generalization of Newtonian gravity in 
terms of an evolution equation for the Eij in addition to the Poisson equation (| ) which 
then acquires a status of a constraint equation. As it has been pointed out in Section ^, 
there are a number of approaches known to modify the Newtonian gravity satisfying the 
Poisson equation, see []TS[, []2D|, |]2T| and references therein. One possibihty is to accept 



the so-called local tidal approximation 

The evolution equations for the tidal force tensor Eij 

dE- ■ 

_^ + OE.j - - 3(5'='afc(,E,), - 5^'ujki^E,)i + A'KGpa,, = (70) 



which is a generalization of the Zel'dovich approximation 39|, |0|| by taking into account 
the second spatial derivatives in the equation of motion for a fluid particle With 



such an evolution equation (^) for the tidal force tensor Eij (^) added, the system of 
the Navier-Stokes-Poisson equations in kinematic quantities (^), (H), (0), (pTD-(^), (|59D- 
(|6T1) , (^)- (|66|) and (|69D becomes closed. Then one can look for a solution of the system 
with a suitable set of initial and boundary conditions characteristic for the problem under 
study. One can assume another evolution equation for the tidal force tensor Eij which 
may be dictated by a particular Newtonian cosmological model. 

(B) The system is overdetermined since it has been derived by differentiation of the 
Navier-Stokes-Poisson equations (^, (||), and (p^Ol). It now contains higher derivatives 
of the fluid velocity through the expansion scalar 9 = 9{x^,t) (1-59), the shear tensor 
aij = (7ij{x'',t) (1-58) and the vorticity vector u;* = Lj\x'^,t) (1-50) and these kinematic 
became now additional unknowns to the velocity field = u^{x\t), the fluid density 
scalar p = p{x^,t) and the gravitational potential scalar = (j){x\t). There are 14 
unknowns, namely, 3 components of the fluid velocity u*, 1 function of the density p, 1 
function of the gravitational potential 0, 1 function of the expansion 6, 5 components of 
the trace-free symmetric shear tensor aij and 3 components of the fluid vorticity vector 
u!\ All together they should satisfy 69 equations, namely, 9 evolution equations (|5TD - 
(|5BD, the equation of continuity (pUj), the equation of state (g), the Poisson equation (|), 
5 evolution equations (0), 9 constraints including 3 constraints (|59|), 1 constraint (|60|), 
5 constraints (pT]), 16 integrability conditions including 3 integrability conditions (0), 5 
integrability conditions (|65| ) and 8 integrability conditions (|66|) and 9 constraints for the 
velocity (^) which satisfy 18 integrability conditions (^ ) and (|5B|). 

(C) Even after the reformulation in terms of the kinematic quantities of the fluid 
expansion scalar 6 = 6{x^,t) (1-59), the fluid shear tensor a^j = (Jij{x^,t) (1-58), and the 
fluid vorticity tensor Uij = uJij{x^,t) (1-49) or the fluid vorticity vector cu* = to'*(x-', t) (I- 
50) the system of the Navier-Stokes-Poisson equations (|), (|), (||), (I5l|)-(|D, (HD-dl]), 

and (|7Dp do contain the fluid velocity u\x^ ,t) explicitly. Indeed, 



the material derivatives in the evolution equations (0)-(0) have the inertial term Ui^u-' 
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involving u'^{x^ ,t). That means from the mathematical point of view that the information 
stored in the spatial derivatives of the fluid velocity Uij{x^,t) through the decomposition 
(|5^ is not enough to fully specify the fluid velocity field u^{x\t) in terms of the kinematic 
quantities. Therefore, if only the kinematic quantities 6{x'^,t), ijij{x^,t) and uj^{x\t) are 
considered to be the unknowns, the system of the Navier-Stokes-Poisson equations is 
actually a system of integro-differential equations. 

(D) The space of solutions of the system (|), (|), ([ID, (|TD-(|6D, (il)-(|6l|), (H"®, 
and (|7y) is bigger than that of the Navier-Stokes-Poisson equations 



(H) and (|ToD . This is reflected in the presence, first of all, of an evolution equation for the 
tidal force tensor (|70|) which does not take place in the original system. Secondly, there 
are constraints and integrability conditions which include even higher derivatives of the 
tensor Wj j. Since the decomposition law (|5^ is a Pfaffian system when it is considered as 
a system of nonhomogeneous partial differential equations for the 3 unknown components 
of velocity Ui, due to the Frobenius theorem the integrability conditions ( ^51) 



and ( pBD are necessary and sufficient for the existence of a solution locally. The problem, 
however, is that the identities (^) and (|56D taken as a system of nonhomogeneous partial 
differential equations for the tensor Uij form a system of partial differential equations 
which is not Pfaffian. In such a case, their integrability conditions (|6^ ) -(l66| ) are only 
sufficient conditions for integrability of (|55D and (|56D in general (see, for instance, |^ and 
references therein). One may expect they become also necessary for some cosmological 
fluid configurations. When and under which assumptions about the physics of a self- 
gravitating fluid, and finally, for which class(es) of Newtonian cosmologies this can occur 
definitely deserves a careful study. It should be pointed out here that this problem of the 
integrability of the Navier-Stokes-Poisson equations in kinematic quantities is affected also 
by the structure of the evolution equations for the tidal force tensor Eij ([70|), or another 
equation of this kind, which has been added to close the system. 

(E) The dynamics of the self-gravitating fiuid governed by the system of Navier-Stokes- 
Poisson equations in terms of kinematic quantities (|p, (^), (|I^), (PD-(pBD, ), 



d6D, (^), p8| ) and (ffOD is not, strictly speaking, due to Newtonian gravity of the 
cosmological fiuid, but rather due to its generalization in accordance with the evolution 
equation ( |70D for the gravitational potential through the tidal force tensor Eij (|28|). The 
universes corresponding to such dynamically evolving cosmological fiuid configurations 
under this modified Newtonian gravitational force will be still called Newtonian, since 
if the Newtonian cosmological principle is satisfied, the tidal force tensor should van- 
ish asymptotically and ( pUD on the largest cosmological scales and the Newtonian 
universes homogeneous and isotropic (^2]) on such large scales are governed then by the 
Newtonian gravitational potential satisfying the Poisson equation (H). If a Newtonian uni- 
verse is inhomogeneous and/or anisotropic the full system of the Navier-Stokes-Poisson 
equations in terms of kinematic quantities must be taken into consideration. 

Though the system of the Navier-Stokes-Poisson equations in terms of kinematic quan- 
tities look a bit complicated as compared with the original Navier-Stokes-Poisson equa- 
tions, the former has proved to be mathematically efficient and physically adequate in 
cosmological studies, in particular, where some simplifying assumptions can be adopted 
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see 



and recent applications |^^-|^ and references therein). 



12 The Raychaudhuri Equation in a Homogeneous 
and Isotropic Newtonian Universe 



The evolution equation for the expansion scalar (|5l|) is known to have the very similar form 
to the Raychaudhuri equation for the expansion scalar of the 4- velocity of the cosmological 
fluid in general relativity with the only different term with pressure involved in 



the general relativistic case. This equation is very significant in both Newtonian and 
general relativistic cosmologies since the expansion scalar of the evolving cosmological 
fluid serves as its fundamental parameter. For an homogeneous and isotropic Newtonian 
universe it is the only dynamical fluid characteristic. One can show in this case that 
the Raychaudhuri equation (^Tj) is equivalent for the Friedmann equation (^), and the 
expansion scalar is then directly expressible in terms of scale factor [O], [|12|, ||37[|. 



Theorem 7 (The Homogeneous and Isotropic Raychaudhuri Equation) The ev- 
olution of a Newtonian universe satisfying the Newtonian cosmological principle is gov- 
erned by the Raychauhuri equation 

^ + -9'' + 4nGp - A = (71) 
at 3 

for the cosmological fluid expansion scalar 9 = 9(t), while the shear tensor aij{x^,t), 
vorticity vector uj\x^ ,t) and the tidal force tensor Eij{x^ ,t) vanish 

= 0, = 0, Eij = 0. (72) 

Then the Raychauhuri equation is equivalent to the Friedmann equation ^3^ ) for the fluid 
scale factor R = R{t) with 

Proof. Let us consider a Newtonian universe flUed with a cosmological fluid satisfy- 
ing the Newtonian cosmological principle and governed by the system of Navier-Stokes- 
Poisson equations in terms of kinematic quantities (|), (|), (0), (|5l|)- ([56|) , (|59D -(|6lD, 



(P^) and (|7Dp. The same analysis as in the case of an homogeneous and 
isotropic Newtonian universe satisfying the Newtonian cosmological principle, Section 
shows that the fluid density p{x^,t) and fluid pressure p{x'',t) are functions of absolute 
time t only and the velocity of the cosmological fluid M*(x^,t) has the following form, 

p = p(t), p = pit), u\t) = H{t)x\ (74) 

Calculation of the tensor Uij and the kinematic quantities for the fluid gives 

= H{t)6ij, 9 = H{t), Oij = 0, J = 0. (75) 
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The kinematic quantities ([7^) together with a sohition to the equation of continuity (^Op 
satisfy the equations (^), (|53| ) -(|56| ) , (|5^ )-(|61|), (|6^) and (|69| ) from the system of Navier- 
Stokes-Poisson equations in terms of kinematic quantities. The evolution equation ( ^21) 
for the shear tensor results in everywhere vanishing tidal force tensor Eij for the fluid 
(^31) which satisfies the evolution equation for the tidal force tensor (^). The vanishing 
tensor Eij ( ^5D is consistent with the solution of the Poisson equation (|]). Then the 
Raychaudhuri evolution equation for the expansion scalar 6 = H{t) is the only remaining 
equation which takes the form (^). By taking into account the definition of the scale 
factor R{t) of the cosmological fluid (^), see also the representation of the expansion 
scalar (1-61), one can easily derive the evolution equation ( ^5|) for the scale factor from 
the Raychaudhuri equation (|7lD . The first integral of ( ^5]) leads the Friedmann equation 
for the scale factor R{t). QED 



This result shows, in particular, that for the case of homogeneous and isotropic New- 
tonian universes the system of Navier-Stokes-Poisson equations in terms of kinematic 
quantities (|), (|), ([T|), (|5l|)-(|5|), (Hl-dlD, (H-®, ®, (H) and ® is consistent. 
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